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Executive  Summary 

As  the  U.S.  Navy  aggressively  pursues  the  design  of  high-speed  aluminum  vessels  for  littoral 
combat  operations,  there  is  a  growing  need  for  on-board  hull  monitoring  systems  capable  of 
detecting  hull  deterioration  and  damage.  In  particular,  fatigue  damage  is  anticipated  due  to  the 
combination  of  non-traditional  hull  forms  and  demanding  seaway  loads.  Fatigue  is  an 
accumulated  degradation  process  that  occurs  under  cyclic  loading,  eventually  inducing  cracking 
at  stress  concentration  points.  In  this  project,  a  strategy  is  proposed  for  fatigue  life  estimation  of 
a  ship  hull  using  a  wireless  sensor  network  installed  in  the  hull  for  autonomous  health 
monitoring.  Specifically,  the  computation  and  storage  capabilities  of  each  wireless  sensing  unit 
are  leveraged  to  process  large,  continuous  amounts  of  raw  strain  gage  data  for  fatigue  life 
information,  thereby  presenting  the  user  with  very  little  data,  yet  very  meaningful  decision 
making  information.  First,  a  rainflow  counting  algorithm  is  implemented  as  a  continuous  time- 
domain  approach  to  fatigue  estimation.  A  peak  picking  algorithm  reduces  time  history  data  to  a 
minimum  set  of  data  points  that  are  needed  for  fatigue  analysis.  Cycles  are  defined  by  a  raintlow 
counting  procedure  and  are  kept  on-board  the  wireless  sensing  unit  in  a  histogram  structure  for 
long-term  storage.  Fatigue  damage  is  accumulated  by  the  Palmgrcn-Miner  linear  summation 
method.  Second,  a  statistical  approach  called  the  Dirlik  procedure  is  used  which  relies  on 
spectral  moments  of  the  stress-time  history  power  spectral  density  (PSD)  function.  The  main 
motivation  for  this  approach  is  to  use  a  few  statistical  parameters  from  the  power  spectrum  to 
estimate  fatigue  life  rather  than  doing  the  continuous  collecting  and  processing  that  is  conducted 
in  the  rainflow  procedure.  The  Dirlik  procedure  outputs  a  probability  density  function  (PDF)  of 
stress  ranges,  tailored  to  resemble  the  raintlow  counting  results.  The  PDF  is  converted  to 
damage  using  an  S-N  curve  and  again  accumulated  by  the  Palmgren-Miner  method.  Damage 
results  arc  kept  on-board  for  long-term  storage  aw  aiting  user  query.  Communication  of  results  is 
done  only  by  request,  dramatically  reducing  the  power  consumption  in  the  wireless  sensing 
network. 

Experimental  tests  are  conducted  on  an  aluminum  hull  stiffened  element  specimen  as  part  of  the 
Monitored  Aluminum  Hull  Integrity  (MAHI)  test  program  to  verify  the  embedded  fatigue  life 
estimation  procedures.  Testing  is  carried  out  at  the  Naval  Surface  Warfare  Centcr-Carderoek 
Division  (NSWCCD)  Fatigue  Testing  Laboratory.  Among  the  key  objectives  of  the  MAHI  test 
program  is  the  demonstration  and  evaluation  of  SUM  technologies.  A  tethered  strain  gage 
monitoring  system  with  16  channels  is  employed  for  laboratory  data  collection.  A  network  of 
Narada  wireless  sensor  units  is  used  to  collect  multiple  strain  gage  channels  performing 
continuous  data  collection  as  well  as  the  embedded  rainflow  and  Dirlik  fatigue  estimation 
procedures.  Strong  agreement  is  found  between  the  Narada  and  NSWCCD  data.  Furthermore, 
the  accuracy  of  the  embedded  fatigue  estimation  algorithms  is  confirmed. 
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1.  Introduction 

In  recent  years,  the  need  for  high-speed  littoral  combat  ships  (LCS)  has  driven  U.S.  Navy  interest 
in  aluminum  hulls  of  unconventional  design  (e.g.,  multi-hull  catamaran  style  ships).  For 
example,  the  Joint  High-Speed  Vehicle  (JHSV)  and  Littoral  Combat  Ship  (LCS)  programs  are 
examples  of  major  ship  design  and  procurement  programs  instituted  by  the  U.S.  Navy  (Hess 
2007).  While  aluminum  allows  the  ship  weight  to  be  dramatically  reduced  as  compared  to  the 
use  of  steel,  aluminum  introduces  a  number  of  significant  challenges  for  long-term  ship 
operation  and  management  (O&M).  Foremost  among  those  concerns  is  fatigue  cracking  in  the 
aluminum  hull.  Fatigue  is  predicted  to  be  a  serious  O&M  issue  due  to  a  general  lack  of  design 
and  analytical  methods  for  aluminum  hulls,  especially  hulls  that  will  experience  the  excessive 
seaway  loads  anticipated  for  LCS  (Tveiten.  et  al.  2007).  Since  the  mid-1950’s,  the  Navy  has 
utilized  aluminum  in  the  design  and  construction  of  the  deckhouses  of  their  steel  destroyer-class 
ships  (Siclski  2007).  However,  fatigue  issues  have  emerged  with  these  deckhouses,  especially 
the  rate  of  growth  of  fatigue  cracks  resulting  in  severe  deckhouse  damage  in  extreme  sea 
conditions.  As  the  Navy  moves  toward  aluminum  for  the  construction  of  LCS  hulls,  it  is 
interested  in  the  adoption  of  long-term  hull  monitoring  systems  that  can  closely  monitor  hull 
responses  to  seaway  loads.  The  benefits  of  hull  monitoring  are  enormous.  First,  hull  response 
data  can  be  used  to  validate  design  assumptions  leading  to  improvements  in  the  design  of  future 
vessels.  Second,  hull  response  data  would  be  integral  to  the  health  assessment  of  the  ship  (i.e., 
structural  health  monitoring  (SHM))  and  guide  O&M  decision  making  of  the  port  engineer.  The 
early  detection  of  hull  damage  and/or  deterioration  allows  the  Navy  to  resolve  structural  issues  in 
a  more  cost-effective  and  timely  manner  leading  to  dramatic  reductions  in  the  total  life-time  cost 
of  ownership  of  a  naval  vessel. 

Researchers  have  begun  to  investigate  the  development  of  structural  health  monitoring  (SHM) 
systems  for  high-performance  ships  (Farrar,  et  al.  2009).  Any  SUM  system  designed  and 
deployed  for  aluminum  hull  ships  must  exhibit  the  ability  to  detect  the  onset  of  fatigue  in  the  hull 
components.  This  project  explores  data-driven  algorithms  that  can  be  autonomously  executed  by 
the  hull  monitoring  system  to  identify  fatigue  damage.  Known  hull  hot-spots  ( e.g .,  base  plates, 
heat  affected  zones  near  welds)  arc  often  instrumented  with  strain  gages  to  record  the  hull 
response  under  cyclic  wave  loading  (Bachman,  et  al.  2007).  The  strain  gages  would  normally  be 
interfaced  to  a  wired  hull  monitoring  system.  However,  in  recent  years,  wireless  telemetry  has 
been  explored  in  the  shipboard  environment  because  of  their  potential  to  eliminate  monitoring 
system  wires  (Lynch,  et  al.  2009).  The  reported  benefits  of  wireless  hull  monitoring  include 
lower  installation  costs  and  the  installation  of  a  larger  number  of  sensors  in  the  hull.  Regardless 
of  the  type  of  hull  monitoring  system  employed  (e.g.,  wired  or  wireless),  a  scalable  approach  to 
processing  hull  response  data  for  an  accurate  assessment  of  the  hull  health  is  still  lacking.  Until 
such  algorithms  are  shown  accurate  and  reliable  in  the  detection  of  hull  deterioration,  acceptance 
of  SHM  in  the  naval  engineering  community  will  be  limited. 


2.  Technical  Objectives 

This  project  explores  data-driven  algorithms  that  can  be  autonomously  executed  by  the  hull 
monitoring  system  to  specifically  identify  fatigue  damage.  The  approach  to  fatigue  analysis 
depends  on  use  of  hull  strain  time  history  data  to  automate  the  process  of  damage  detection. 
Time-  and  frequency-domain  algorithms  arc  adopted  to  continuously  and  autonomously  count 
strain  peaks  in  the  recorded  strain  time  histories.  Using  peak  amplitude  strains  (S)  and  the 
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Figure  1 :  Rainflow  cycle  analysis  via  closed  loop  identification 


accumulated  number  of  peaks  (N),  the  Palmgren-Miner  rule  and  the  appropriate  S-N  curve  for 
the  ship’s  aluminum  alloy  arc  selected  to  roughly  track  the  fatigue  life  of  the  instrumented  hull 
detail.  The  selection  of  these  methods  of  tracking  the  fatigue  life  of  hull  components  is 
motivated  by  their  computational  simplicity.  In  particular,  these  fatigue  estimation  algorithms 
are  embedded  into  the  computational  cores  of  Narada  wireless  sensor  nodes;  Narada  has  been 
previously  validated  for  hull  monitoring  on  the  FSF-1  ScaFighter  (Lynch,  et  al.  2009). 

The  specific  project  objectives  are: 

•  Rain  flow  assessment  of  strain  time  histories  will  be  implemented  along  with  Miner’s 
rule  in  the  computational  core  of  a  wireless  sensor  for  fatigue  monitoring  of  aluminum 
components.  Experimental  validation  will  be  performed  on  aluminum  coupon  specimens 
loaded  under  uniaxial  load. 

•  Spectral  methods  to  peak  extraction  will  also  be  implemented  with  Miner’s  rule  in  the 
computational  cores  of  a  wireless  sensor.  Again,  experimental  validation  will  be 
conducted  using  aluminum  coupons  loaded  under  uniaxial  load. 

•  Third,  the  embedded  fatigue  analysis  will  be  validated  experimentally  on  complex 
aluminum  hull  connections  at  NSWC-Carderock. 


3.  Theoretical  Background  and  Algorithms 

3. 1  Rain  Flow  Peak  Extraction  and  Miner 's  Rule  for  Fatigue  Estimation 

Fatigue  life  monitoring  begins  by  identifying  cycles  within  a  complex  load  time  history.  A 
number  of  cycle  counting  methods  are  used  to  reduce  irregular  load  histories  into  a  collection  of 
constant  amplitude  events  such  as  rainflow  counting  (primarily  used  herein),  range-pair  counting, 
and  racetrack  counting.  Rainflow  cycle  counting  has  shown  to  be  among  the  superior  methods 
for  cycle  counting  of  irregular  loads.  In  the  rainflow  method,  cycles  are  identified  in  a  manner  in 
which  closed  hysteresis  loops  arc  identified  from  the  stress-strain  response  of  a  material  subject 
to  cyclic  loading.  As  shown  in  Fig.  1,  closed  hysteresis  loops  can  be  identified  from  the  Strain- 
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Figure  2:  Histogram  of  cycle  accumulation  Figure  3:  Strain-life  curve 

time  history  shown.  Ranges  A-D,  B-C,  E-F,  and  G-H,  would  represent  cycles  counted  under 
rainflow  counting  techniques.  Rainflow  counting,  however,  is  originally  intended  to  be  carried 
out  once  the  entire  strain  history  is  known,  since  counting  starts  and  ends  at  the  maximum  peak 
or  valley.  Due  to  the  limited  memory  available  in  wireless  sensors,  it  is  not  possible  to  wait  until 
the  entire  load  history  has  been  realized  before  fatigue  accumulation  can  be  calculated.  Rather, 
the  ‘one-pass’  rainflow  counting  algorithm  can  be  alternatively  used  since  it  does  not  require  the 
full  time  history  record.  This  method,  which  was  used  for  embedment  in  this  study,  is  a  vector- 
based  counting  algorithm  first  demonstrated  by  Downing,  et  al.  (1976)  and  later  modified  by 
Okamura,  et  al.  (1979)  to  account  for  half  cycles.  The  embedded  rainflow  counting  algorithm 
starts  by  arranging  sampled  strain  time  history  measurements  into  a  single  vector.  From  the  set 
of  strain  data,  peak  and  valley  points  inherent  in  cyclic  measurements  arc  identified  and  stored. 
The  typical  rainflow  counting  procedure  is  then  performed  on  the  set  of  peaks  and  valleys, 
identifying  closed  hysteresis  loops  and  logging  those  ranges  as  cycles.  For  each  cycle,  both  the 
strain  amplitude  and  mean  strain  arc  recorded.  This  rainflow'  counting  procedure  identifies  the 
same  cycles  as  the  traditional  rainflow  procedure  that  uses  the  entire  strain  time  history.  This 
fact  makes  the  ‘one-pass’  rainflow'  counting  procedure  very  attractive  for  continuous  real-time 
monitoring  of  fatigue  life  using  wireless  sensors. 

The  ‘one-pass’  rainflow'  counting  technique  can  be  operated  in  a  real-time  manner  for  continuous 
monitoring  of  fatigue  life.  The  compression  of  a  strain  time  history  from  sampled  points  to 
cycles  though,  is  not  enough  to  monitor  fatigue  in  the  long-term.  Further  data  compression  is 
realized  by  accumulating  identified  cycles  of  mean  and  amplitude  into  a  histogram  similar  to  the 
one  shown  in  Fig.  2.  A  fixed  size  histogram  allows  for  a  priori  allocation  of  the  available 
memory  integrated  w'ith  the  wireless  sensor,  allowing  for  the  near  perpetual  accumulation  of 
fatigue  cycles  and  thus  continuous  “real-time”  fatigue  life  monitoring.  Considerations  of  bin 
sizing  in  the  mean  and  amplitude  axes  will  be  of  importance  in  terms  of  mean  and  amplitude 
accuracy,  as  well  as  in  terms  of  the  scarce  memory  consumed. 

Strain-life  methods  which  account  for  mean  strain  effects  arc  used  in  this  project  to  predict 
fatigue  life.  First,  the  total  strain  amplitude  can  be  expressed  as  the  sum  of  an  elastic  strain 
amplitude  and  a  plastic  strain  amplitude,  each  represented  linearly  against  cycles  to  failure  on  a 
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log-log  plot.  Total  strain  amplitude  versus  eyelcs  to  failure,  Nt,  before  mean  stress  (or  mean 
strain)  correction  is  written  as: 

At  <*f  „  (1) 

—  =  -f{2N,)h+  ef(2NfY  Kl> 

where:  y  =  strain  amplitude 

<jf'  =  fatigue  strength  coefficient 
E  =  modulus  of  elasticity 
2Nf  =  reversals  to  failure 
e/  =  fatigue  ductility  coefficient 
b  =  fatigue  strength  exponent 
c  =  fatigue  ductility  exponent 

The  constants  b  and  c  represent  the  slopes  of  the  clastic  and  plastic  strain.  Similarly,  <'-L  and  e/ 
are  the  y-axis  intercepts  of  the  elastic  and  plastic  strain  curves  in  Fig.  3.  Eq.  (1)  is  applicable  if 
dealing  with  full  cycles  with  zero  mean  strain.  However,  several  researchers  have  proposed 
modifications  to  the  strain-life  relationship  of  Eq.  (1)  to  account  for  mean  stress  effects  including 
Morrow  (1968)  and  Manson  and  Halford  (1981).  These  relations  all  have  their  own  advantages, 
but  require  the  mean  stress  in  order  to  evaluate  the  number  of  cycles  to  failure.  A  modification 
to  the  strain-life  relation  using  mean  strain  instead  of  mean  stress  would  be  better  because  strain 
is  measured  and  not  stress.  The  need  to  track  stress-strain  around  a  hysteresis  curve  to  determine 
mean  stress  is  taken  out  of  the  picture,  making  the  embedded  process  easier  to  implement.  One 
sueh  empirical  relationship  between  total  strain  range,  Ae,  and  cycles  to  failure  is  adopted  here. 
This  empirical  relationship  which  accounts  for  mean  strain  instead  of  mean  stress  is  w  ritten  as: 

2(1 (2) 
[(4fy-l)(l-R)a  +  (2)a]1/fl 


where.  R  tminf6max  it  |£ma*l  ^  l^mml 

R  tmax/tmin  h  l^maxl  ^ 

—  tmax  —  Emin 


Here,  a  is  equal  to  -1/c.  The  fatigue  duetility  coefficient,  e/,  fatigue  strength  exponent,  b,  and 
fatigue  duetility  exponent,  c,  determined  for  the  strain-life  relationship  are  used  both  directly  and 
in  computing  the  material  constant,  a.  These  parameters  are  best  evaluated  by  performing  eyelie 
testing  in  the  laboratory.  When  fatigue  data  is  not  available  or  easily  obtained,  these  parameters 
can  be  estimated  from  statie  properties  or  by  using  other  estimation  methods. 

In  this  study,  the  empirical  fatigue  law  of  Eq.  (2)  will  be  used.  To  estimate  the  model  parameters 
(fatigue  duetility  coefficient,  e /,  fatigue  duetility  exponent,  c,  and  fatigue  strength  exponent,  b), 
the  uniform  material  law  by  Baumel  and  Seeger  (1990)  is  used.  For  aluminum  and  titanium 
alloys,  the  strain-life  equation  is  estimated  as 

y  =  1.67  y  (2N,)-°9S  +  0.35  (2Nf)-°69  ^ 


Page  8  of  29 


i 


where:  ab  =  yield  stress. 

While  Eq.  (3)  will  not  be  used,  one  can  however  extract  the  fatigue  ductility  coefficient,  e/, 
fatigue  ductility  exponent,  c,  and  fatigue  strength  exponent,  b  from  it.  With  knowledge  of  these 
constants,  the  material  constant,  a,  can  be  estimated  and  Eq.  (2)  can  be  written  in  its  final  form  as: 

= _ 2U-J003S _  (4) 

[(4 Nr  -  1)(1  -  R)~0A  +  (2)_04]1/-0'4 

In  this  form,  the  mean  strains  and  strain  amplitudes  obtained  through  rainflow  counting  arc 
sufficient  for  determining  the  fatigue  life,  Nf,  of  an  instrumented  component.  Although  it  is 
possible  to  embed  a  mean  stress  procedure,  it  will  require  more  work  for  the  wireless  sensor  and 
may  be  unnecessary  as  the  prediction  proposed  in  Eq.  (4)  has  been  shown  to  compare  extremely 
well  with  data  received  from  alloys  tested  under  a  variety  of  different  tensile  and  compressive 
mean  strains  (Collins  1981). 

The  fatigue  life  corresponding  to  each  cycle  is  an  indication  of  the  amount  of  damage  imposed 
on  the  material  due  to  that  specific  cycle.  In  this  way,  one  can  start  to  accumulate  and  monitor 
damage.  The  damage  summation  of  each  cycle  is  done  using  the  Palmgrcn-Miner  linear  damage 
hypothesis  originally  proposed  by  Palmgrcn  (1924)  and  later  modified  by  Miner  (1945).  The 
Palmgren-Mincr  rule  is  written  as: 


k 


(5) 


where:  D  =  accumulated  damage 

k  =  total  number  of  cycles  in  a  loading  spectrum 
i  =  ith  applied  stress/strain  level 
nt  =  number  of  cycles  at  stress/strain  level  i 
Nt  =  fatigue  life  at  stress/strain  level  i 

This  method  implies  that  failure  occurs  when  the  summation  of  cycle  ratios,  n/Nj,  is  equal  to  1 . 
It  should  be  noted  that  by  using  strain-life  methods,  the  initial  cracking  is  predicted  instead  of 
complete  failure.  Further  monitoring  of  fatigue  damage  (/.<?.,  after  the  initiation  of  cracking) 
requires  eraek  propagation  methods  which  are  beyond  the  scope  of  this  project. 

The  value  residing  in  each  bin  of  the  accumulated  cycle  histogram  represents  the  number  of 
cycles  at  a  specific  strain  amplitude  and  mean.  The  position  of  each  bin  determines  the  fatigue 
life  of  that  particular  bin,  since  it  represents  the  strain  amplitude  and  mean  strain  required  for  the 
strain-life  relation.  In  the  embedded  implementation,  cycles  incurring  damage  below  a  specified 
threshold  result  in  their  elimination  from  the  analysis.  Signal  noise  will  generate  high  amounts 
of  low  amplitude  cycles,  which  when  accumulated  over  an  excessive  amount  of  time  may  falsely 
contribute  to  the  expended  life  of  the  material.  Although  certain  materials  have  no  defined 
fatigue  limit,  this  liberty  is  assumed  safe  as  noise  level  amplitudes  are  far  below  reasonably 
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assumed  estimations  for  fatigue  limits  of  these  materials.  It  is  important  to  use  damage  as  the 
parameter  for  which  eyeles  arc  eountcd  or  eliminated  from  analysis.  Very  small  amplitude 
eyeles  at  very  high  mean  stress  may  result  in  relevant  fatigue  damage,  and  should  not  be 
eliminated  from  analysis. 

3.2  Spectral  Methods  for  Fatigue  Estimation 

A  frequency  domain  approach  to  eyele  identification  for  estimation  of  fatigue  accumulation  is 
also  explored  as  an  alternative  to  the  proposed  rainflow  methods.  Sinee  many  loading 
environments  arise  as  stochastic  processes  (t?.g.,  live  loads,  wind,  waves),  it  is  appropriate  to  use 
power  speetral  density  (PSD)  as  a  statistical  function  describing  operational  loads.  The  main 
motivation  for  the  frequency  domain  approach  is  to  use  only  a  few  statistical  parameters  dedueed 
from  a  power  speetrum  to  estimate  fatigue  life  rather  than  doing  continuous  eyele  counting  in  the 
time  domain  as  is  the  ease  in  rainflow  counting  methods.  Many  results  exist  for  the  ease  of 
narrow  band  stationary  and  ergodie  random  signals,  especially  those  having  a  Rayleigh  stress 
range  distribution;  however,  one  would  like  to  have  a  procedure  capable  of  handling  wide  band 
loading  as  well  (Kirkby  1966).  The  Dirlik  (1985)  approach  used  here  provides  a  closed  form 
expression  linking  any  given  power  speetral  density  to  a  probability  density  function  of  rainflow 
ranges  when  the  random  process  is  stationary  and  ergodie.  The  Dirlik  procedure  uses  a  Monte 
Carlo  approach  to  generate  stress  time  histories  from  various  PSD  shapes.  Time  domain 
rainflow  counting  is  performed  and  a  probability  density  function  of  rainflow  ranges  is  linked  to 
the  corresponding  PSD  shape.  Dirlik’s  solution  is  empirically  derived  by  fitting  rainflow  range 
distributions  to  simulated  data  and  is  technically  restricted  to  Gaussian  processes.  Theoretically 
based  solutions  exist  under  eertain  conditions,  (e.g..  Bishop  and  Sherrat  (1990))  but  Dirlik’s 
method  remains  among  the  most  frequently  employed  techniques  due  to  its  simplicity  and  due  to 
its  sufficient  level  of  accuracy  under  required  conditions. 

The  bulk  of  signal  processing  in  Dirlik’s  method  lies  in  the  Fourier  transform  used  to  convert 
time-history  strain  data  to  the  frequency-domain.  A  fast  Fourier  transform  (FFT)  is  used  to 
determine  the  power  spectral  density  (PSD)  function  of  strain  data  collected.  The  PSD  is  the 
base  of  all  calculations  in  the  analytical  procedure.  Specifically  of  interest  are  the  speetral 


moments,  where  the  n'h  moment  is  defined  as; 


co 


(6) 


or,  in  the  ease  of  a  discrete  PSD,  can  be  visually  observed  in  Fig.4  and  written  as 


(7) 


^  fnG(f)df 


Many  probabilistic  eharaeteristies  of  the  signal  arise  from  these  moments.  In  1945,  Riee  (1945) 
provided  the  relationship  for  the  number  of  upward  zero  crossings  and  peaks  per  second  based 
solely  on  speetral  moments,  paving  the  way  for  a  fatigue  life  estimation  based  on  power  speetral 
density.  The  relations  are  seen  below  as: 
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fmox 


Expected  number  of  zero  crossings  per  unit  time,  A0: 


Expeeted  number  of  peaks  per  unit  time,  jU: 


^  = 


Irregularity  faetor,  y: 


Aq  _  m22 

M  Jmnmi 


(8) 


(9) 


(10) 


Irregularity  factor  is  the  measure  of  the  expeeted  ratio  of  zero-crossings  to  peaks,  approaching 
zero  for  the  wide  band  case  and  one  for  the  narrow  band  case.  Other  important  characteristics 
arising  from  the  spectral  moments  include  the  root  mean  square  (RMS),  ax: 


ax  = 


Jmo 


and  mean  frequency,  Xm : 


m, 

Xm  =  — 

mo . 


m2 

m4 


on 


(12) 
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From  the  above  variables,  Dirlik’s  rainflow  -  range  probability  density  function,  Prr,  is  given  as: 


where: 


-Zj  Z  -z2/  ,  -z2/ 

?  /f+C2^ 2e  2a  +  C^ze  f2 

(13) 

X 

Z  =  — 

2  Ox 

(14) 

2 (xm  -  r2) 

Cl  i  +  r2 

(15) 

V-Xm-Cj2 

i  -  r-ci  +  ci2 

(16) 

i  -  r  -  ci  +  c2 
l~  1  —  a 

(17) 

C3  =  1  -  Cl  -  c2 

(18) 

1.25  (y-c3-c2) 

(19) 

The  probability  density  function  of  rainflow  ranges  relates  fatigue  life  to  a  measured  power 
spectral  density.  Expected  damage  can  be  calculated  for  any  kind  of  S-N  relation  by  the 
expected  damage  per  unit  time,  H( D): 


E(D)=ftf 


N(S) 


PRR(S)dS 


(20) 


where  N(S)  is  the  number  of  cycles  to  failure  for  a  stress  range  S.  In  the  case  of  a  straight  line 
log-log  scale  representation  of  the  S-N  curve,  the  following  can  be  used: 


(21) 


where  K  and  b  arc  material  constants.  Substituting  Eq.  (21)  into  Eq.  (20)  for  the  total  damage 
expected  from  time  t= 0  to  T,  the  following  is  obtained: 


oo 

E(D)t=TK-'ii  j  ShPRR(S)dS 

0 


(22) 
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Figure  5:  Narada  wireless  sensing  node:  (top)  key  components  highlighted; 
(bottom)  small  network  of  nodes  ready  for  deployment. 


4.  Experimental  Validation  of  the  Embedded  Fatigue  Estimation  Algorithms 

4. 1  Implementation  in  the  Narada  Wireless  Sensor  Node: 

The  aforementioned  fatigue  estimation  algorithms  were  selected  for  embedment  in  the 
computational  core  of  the  Narada  wireless  sensor  node.  This  particular  wireless  sensor  node  has 
been  previously  validated  during  sea  trials  of  the  FSF-I  SeaFighter  LCS  (Lynch,  et  al.  2009). 
The  Narada  wireless  sensing  unit  shown  in  Fig.  5  was  developed  at  the  University  of  Michigan. 
The  Narada  uses  an  Atmel  Atmegal28  microprocessor  with  128kB  of  external  SRAM  for  data 
storage  and  computation.  The  external  memory  allows  the  unit  to  store  up  to  64,000  data  points 
at  one  time.  Wireless  communication  is  realized  via  a  Chipeon  CC2420  IEEE802.15.4- 
eompliant  wireless  radio,  making  the  unit  exceedingly  versatile  for  developing  large,  sealable 
wireless  sensor  networks.  This  unit  utilizes  a  four  channel,  16-bit  Texas  Instalments  ADS8341 
ADC  for  data  acquisition,  and  a  two  channel,  12-bit  Texas  Instruments  DAC7612  digital-to- 
analog  converter  (DAC)  for  actuation  capability. 

The  embedded  rainflow  fatigue  estimation  procedure  utilizes  the  128kB  of  external  SRAM  on 
the  Narada  wireless  sensor  node.  The  external  memory  bank  is  divided  into  two  64kB  halves 
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Figure  6:  Memory  map  for  rainflow-based 
fatigue  monitoring  procedure 


Figure  7:  Memory  map  for  embedded 
Dirlik  fatigue  monitoring  procedure 


denoted  as  low  and  high  (Fig.  6).  In  the  rain  flow  implementation,  strain  data  acquired  is  stored 
on  the  low  side  in  three  separate  blocks  with  a  eapaeity  of  nearly  20kB  eaeh.  This  allows  for  the 
acquisition  and  processing  of  data  sets  containing  nearly  10,000  data  points  of  2-bytes  (or  16- 
bits).  The  three  stage  data  acquisition  guarantees  the  continuous  acquisition  of  strain  data  in  one 
block  while  fatigue  life  procedures  are  performed  simultaneously  on  a  previously  attained  set  of 
strain  data.  After  a  stack  of  memory'  is  filled,  a  peak  picking  algorithm  is  performed  on  the  set  of 
strain  data.  A  set  of  only  peaks  and  valleys  are  stored  in  a  speeifie  location  on  the  high  side  of 
the  external  SRAM.  The  extreme  condition  considered  for  sizing  memory  allocation  occurs 
when  all  points  acquired  during  sampling  result  in  peaks  or  valleys  of  the  strain  response  signal, 
requiring  an  equally  sized  20kB  memory  allocation.  The  rainflow  counting  procedure  is  carried 
out  on  the  set  of  peaks  and  valleys  and  stored  on  the  high  side  of  the  external  SRAM.  Extreme 
conditions  here  would  sec  less  than  half  the  number  of  peak/valley  points  as  the  maximum 
number  of  eyeles,  requiring  lOkB  of  memory  allocation  for  strain  amplitude  and  mean  strain 
eaeh.  The  remaining  20kB  of  external  SRAM  on  the  high  side  is  used  as  a  continuously  updated 
histogram.  In  the  ease  where  the  number  of  eyeles  is  represented  as  a  2-byte  integer,  10,000 
unique  strain  amplitude-mean  strain  combinations  are  available.  Any  m-by-n  product  of 
amplitude  bins,  m,  by  mean  bins,  /?,  less  than  10,000  are  admissible.  In  this  projeet,  eyeles  are 
accumulated  in  a  2-byte  memory  slot,  limiting  the  maximum  number  of  cycles  that  can  be 
accumulated  for  a  speeifie  strain  amplitude  and  mean  strain  to  65,535  eyeles.  All  processing 
procedures  carried  out  on  the  high  side  of  the  external  SRAM  allow  for  continuous  strain  gage 
data  acquisition  on  the  low  side  by  way  of  an  interrupt  on  the  analog-to-digital  converter.  The 
external  SRAM  layout  for  the  rainflow  method  is  summarized  in  Fig.  6. 


The  embedded  spectral-based  fatigue  estimation  is  also  embedded  into  the  j Narada  wireless  node 
external  memory  (Fig.  7).  Again,  the  memory  map  is  laid  out  into  two  64kB  halves.  A 
maximum  of  32,000  2-byte  strain  gage  samples  can  be  collected  at  one  time  but  the  FFT 
computation  restricts  us  to  4096  samples  as  implemented  on  Narada :  henee,  8.2  kB  is  allocated 
for  sampled  strain  data  storage.  Once  the  FFT  algorithm  is  run  with  the  sampled  data,  a  complex 
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Table  1:  Analysis  of  communication  requirements  of  centralized  and  proposed  decentralized 

fatigue  life  monitoring  methods 


Methods 

Transmission  payload  byte 

Centralized 

R  a  Inflow'  counting  performed  on  central  server  after  all  time 
history  data  is  received 

Ndata  *2  byte 

=  30000  *  2  =  60  kbyte 

Decentralized 

Ra inflow'  counting  conducted  on  w  ireless  sensing  nodes 
with  cycle  histogram  sent  to  sener 

Histo  *  2  byte 

=  10000  *  2  =  20  kbyte 

Transmission  reduction  =  ~  67  % 

Note:  Time  history  data  length,  Ndata  =  30000  points 


Histogram  size,  Histo  =  10000  bins 

floating  point  number  is  stored  for  each  point  resulting  in  33  kB  of  FFT  data.  The  subsequent 
PSD  calculation  requires  half  the  data  sizes  as  the  floating  point  numbers  are  real  only.  The 
stress  range  PDF  is  given  a  size  of  2,000  points,  ranging  from  0  to  20  ksi.  It  is  up  to  the  user  and 
available  memory  as  to  how  resolute  the  stress  range  PDF  should  be.  The  sampled  strain  data 
and  corresponding  PSD  are  stored  in  the  lower  half  of  memory  while  the  stress  range  PDF  and 
strain  FFT  are  stored  in  the  higher  half  as  shown  in  Fig.  7. 

By  embedding  the  fatigue  life  monitoring  process  into  Narada ,  great  savings  in  communication 
can  be  realized.  Table  1  shows  an  analysis  of  communication  requirements  of  a  centralized 
counting  implementation  and  a  decentralized  computing  implementation.  The  example  in  Table 
1  is  consistent  with  10  minutes  of  continuous  data  acquisition  at  a  sampling  frequency  of  50  Hz. 
The  decentralized  implementation  sends  one  histogram  at  the  end  of  the  ten  minutes.  A 
transmission  reduction  of  approximately  67%  is  the  result  of  a  ten  minute  experiment.  Fatigue 
life  monitoring  though,  will  require  much  longer  periods  of  data  acquisition  and  will  further 
exploit  this  transmission  reduction. 

Fig.  8  shows  the  transmission  payload  over  a  5  hour  period  for  both  the  centralized 


centralized  computing  communication  costs 
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Figure  9:  Sample  random  amplitude  load  Figure  10:  Dual  data  acquisition  of  strain  gage 


implementation  (where  raw  data  is  continuously  streamed  over  the  w'irclcss  channel)  and  the 
decentralized  implementation  (where  fatigue  histograms  are  loeally  updated  and  occasionally 
transmitted).  Raw  strain  gage  data  increases  as  a  linear  function  of  time  for  any  given  sampling 
frequency  in  the  centralized  implementation.  The  number  of  unique  strain  amplitude  and  mean 
strain  cycles  sent  by  the  histogram  in  the  decentralized  implementation  remain  fixed,  and  thus 
increase  by  the  same  amount  at  eaeh  request.  Histograms  are  sent  onec  eaeh  hour.  At  the  end  of 
the  5  hour  example,  the  decentralized  computation  produces  1700  less  kilobytes  for  transmission 
than  the  centralized  implementation,  resulting  in  a  94%  transmission  reduction.  Further 
increases  in  time  between  histogram  transmissions  will  result  in  even  greater  transmission 
savings.  By  receiving  compact  amounts  of  meaningful  data  transmitted  upon  request  as  opposed 
to  receiving  vast  amounts  of  unprocessed  data  transmitted  frequently,  faster  and  more  efficient 
analysis  of  the  structural  member  under  observ  ation  can  be  conducted. 

4.2  Aluminum  Bar  Specimens  under  Uniaxial  Tensile  Loading: 

For  this  study,  a  simple  specimen  constructed  of  6061-T6  aluminum  alloy  is  used.  The 
specimens  are  1/8  in  thick  bars  roughly  1 7  in  long  and  3  in  wide.  A  10  mm  (0.4  in)  strain  gage  is 
used  (Tokyo  Sokki  Kenkyujo  Co.,  Ltd.  TML  PFL-10-1 1-1L  120Q)  at  the  mid-seetion  of  the  bar. 
Bonding  areas  arc  removed  of  grease  and  oils  and  lightly  polished  before  adhering  strain  gages. 
Fatigue  loading  is  applied  using  an  MTS  Series  318  electrohydraulic  closed  loop  load  unit.  The 
load  unit  consists  of  two  smooth  vertical  columns  that  are  joined  by  two  stiff  cross  members,  one 
being  fixed  with  an  integrally  mounted  hydraulie  aetuator.  Hydraulic  grips  hold  the  specimen  in 
plaee  during  loading.  A  representative  loading  input  is  shown  in  Fig.  9  (this  is  a  random 
amplitude  loading). 

In  an  effort  to  compare  a  wired  fatigue  analysis  w’ith  the  embedded  analysis  executed  by  the 
Narada  wireless  sensing  node,  the  specimen  strain  gage  is  split  with  its  output  interfaced  to  both 
a  wired  and  wireless  systems  as  shown  in  Fig.  10.  The  strain  gage  lead  wires  arc  first  attached  to 
a  Wheatstone  bridge  eircuit  to  convert  the  gage  resistance  change  into  a  voltage  signal.  The 
voltage  output  of  the  bridge  eircuit  is  first  passed  to  the  input  of  an  op-amp  amplification  circuit 
before  being  passed  on  to  the  Narada  wireless  sensor  node.  The  amplifier  circuit  amplifies  the 
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Figure  1 1 :  (a)  Reordered  amplitudes  over  entire  time  history  (b)  wired  and  wireless  amplitude  overlay 
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Figure  12:  Close  up  view  of  relevant  wired,  wireless,  and  histogram  cycle  outputs 

bridge  voltage  by  a  factor  of  50.  The  output  of  the  bridge  eireuit  is  also  passed  to  a  National 
Instruments  BNC-2110  data  acquisition  system  interfaced  to  a  personal  computer  running 
LabView.  The  sampling  rate  used  by  the  Narada  mode  is  50  Hz  while  that  used  by  the  National 
Instruments  system  is  100  Hz. 

First,  the  embedded  rainflow  counting  method  is  tested  and  validated.  Specimens  are  subjeet  to 
variable  amplitude  uniaxial  tensile  eyelie  loading.  The  loading  pattern  consists  of  100  randomly 
generated  peaks  and  valleys  as  shown  in  Fig.  9.  Load  sets  are  looped  continuously  during  testing 
until  the  end  of  the  test.  To  make  a  comparison  between  the  wired  and  wireless  fatigue 
monitoring  system,  the  wireless  sensing  system  is  set  to  transmit  all  cycles,  giving  strain 
amplitude  and  mean  strain  for  each  eyele  at  the  end  of  each  sampling  bloek.  Stress  peak  and 
mean  histograms  are  transmitted  periodically  during  testing  upon  user  request.  Fig.  11a 
compares  wired  and  wireless  rainflow  eyele  counting  amplitudes,  plotting  all  amplitudes 
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Oamage  Accumulation 


Table  2:  Fatigue  Life  Results 


Total  Cycles 

Cumulative  Damage 

%  difference 

Wired 

9642 

.00273 

Wireless 

9570 

.00256 

6.22 

Binned  Wireless 

9684 

.00259 

5.13 

calculated  during  rainflow  counting,  reordered  lowest  to  highest.  It  can  be  seen  that  the  wired 
system  produces  more  low  amplitude  cycles  than  the  wireless  system.  This  can  be  explained  as 
an  effect  of  splitting  the  strain  gage  to  the  two  independent  DAQ  systems,  which  introduce  an 
increase  in  signal  noise  to  the  wired  signal.  Additionally,  the  wired  DAQ  system  samples  twice 
as  fast  as  the  wireless  system  resulting  in  more  peaks  and  valleys  to  consider  during  rainflow 
counting.  Figure  1  lb  shows  the  amplitudes  read  from  each  system  where  the  signal  noise  is 
removed  from  the  wired  strain  data,  resulting  in  a  more  accurate  comparison  of  the  two  systems; 
similar  amplitude  outputs  are  extracted  by  both  systems.  Fig.  12  provides  a  eloser  look  at  the 
wired  and  wireless  amplitudes  that  overlay  in  Fig.  lib.  Fig.  12  also  includes  the  results  recorded 
in  the  histogram  during  testing.  For  the  division  of  the  histogram  bins,  327  unique  amplitudes 
and  2 1  unique  mean  values  arc  selected.  Damage  accumulation  for  the  three  sets  of  results  is 
shown  in  Fig.  13.  Damage  value  represents  the  damage  accumulated  as  calculated  by  the 
Palmgren-Miner  method  (namely,  Eq.  (5)).  The  maximum  damage  results  in  Fig.  13  are  also 
listed  in  Table  2.  Excellent  agreement  between  the  wired  and  wireless  systems  is  found. 

The  Dirlik-based  spectral  fatigue  damage  estimation  was  also  embedded  in  the  Narada  wireless 
sensor  and  validated.  First,  constant  amplitude  tests  were  run  at  various  stress  ranges  in  an  effort 
to  observe  expected  probability  distributions.  Probability  density  functions  resulting  from  the 
embedded  Dirlik  procedure  for  constant  amplitude  testing  at  ranges  1  -8  ksi  are  shown  in  Fig.  14. 
These  tests  arc  run  to  achieve  a  simple  visual  confirmation  of  the  expected  change  in  the 
distribution  for  various  stress  ranges.  As  would  be  expected,  the  area  under  the  curve  increases 
near  the  actual  stress  range  for  each  individual  test.  To  compare  the  embedded  rainflow  counting 
and  embedded  Dirlik  procedures,  the  samples  obtained  from  strain  gage  response  of  a  laboratory 
fatigue  specimen  cycled  under  narrow  band  loading  with  an  RMS  stress  range  of  6  ksi  is  used. 
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Figure  14:  Probability  distribution  for  varying  stress  ranges  for  an  aluminum  coupon 

Table  3:  Time  and  frequency 

domain  number  of  cycles 

and  expected  damage 

Rainflow 

Dirlik 

N 

D 

E[N) 

E[D] 

Sample  1 

5226 

2.0754e-6 

13852 

1.4256e-6 

Sample  2 

5033 

2.5871  e-6 

13485 

1 .0990e-6 

Sample  3 

5307 

1  0109e-6 

15649 

1 .3732e-6 

Sample  4 

4970 

2.0767e-6 

13672 

1.8415e-6 

Sample  5 

4945 

2.5781  e-6 

14437 

7.7836e-7 

Fig.  15  shows  the  (a)  windowed  time  domain  response  and  (b)  PSD  for  a  100.000  point  time 
history  sampled  at  100  Hz.  Also  shown  are  the  corresponding  stress  range  distributions  via  (c) 
rainflow  counting  in  the  time  domain  and  (d)  Dirlik  procedure  using  spectral  moments.  The 
number  of  total  cycles  from  the  rainflow  counting  method  and  Dirlik  method  are  shown  in  Table 
3  along  with  the  damage  estimates  from  each  method.  The  total  number  of  cycles  by  rainflow 
counting  is  simply  given  by  summing  the  number  of  cycles  in  each  bin  of  the  cycle  histogram. 
For  the  spectral  procedure,  the  expected  number  of  cycles  for  a  time  T  is  given  by: 

E[N]  =  E[P]  ■  T  (23) 

where  E[P]  —  n  =  expected  number  of  peaks  per  second. 
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Figure  15:  (a)  Stress  -  time  history  of  fatigue  loading;  (b)  one-sided  PSD;  (c)  rainflow 
counting  algorithm;  (d)  Dirlik  PDF  of  rainflow  ranges. 
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4.3  Validation  on  the  Monitored  Aluminum  Hull  Integrity  (MAHI)  Phase  I  Specimen: 

The  NSWCCD  Monitored  Aluminum  Hull  Integrity  (MAHI)  program  presents  an  opportunity  to 
further  validate  the  embedded  fatigue  life  estimation  algorithms.  The  MAHI  fatigue  test 
program  is  specifically  focused  on  the  design,  fabrication  and  fatigue  testing  of  aluminum  hull 
stiffened-plate  structural  specimens.  In  this  project,  the  first  phase  specimen  is  tested  which 
consists  of  a  single  stiffened  plate  clement  with  an  intersecting  stiffened  bulkhead  element.  The 
MAHI  fatigue  test  specimen  is  shown  in  Fig.  16.  The  specimen  consists  of  a  5083-H116 
aluminum  plate,  a  tee  stiffener  maehined  from  an  A16061-T6  aluminum  wide-flange  beam 
extrusion,  and  a  bulkhead  detail  made  from  A15083-H1 16  plate,  as  well  as  a  vertical  tec  stiffener 
on  the  bulkhead  maehined  from  the  above  6061 -T6  stoek  extrusion.  The  specimens  were 
fabricated  according  to  generally  aeeepted  shipyard  praetiees  so  as  to  respond  with  typieal  eraek 
initiation  behavior. 

Fatigue  testing  was  performed  in  the  NSWCCD  Fatigue  Testing  Laboratory  in  Carderock,  MD. 
The  NSWCCD  fatigue  laboratory  contains  six  uniaxial  scrvo-hydraulieally  controlled  load 
frames  of  varying  capacity,  all  having  full  tension  and  compression  eyelie  loading  capability.  An 
Instron  220-kip  machine  was  used  in  the  experimental  tests  for  this  study.  The  Instron  fatigue 
testing  maehinc  is  computer-interfaeed  with  analog  controllers  to  allow'  for  complex  test  loadings. 
The  test  machines  arc  also  equipped  with  wired  data  acquisition  system  with  channels  dedicated 
to  recording  applied  load,  displacement,  and  strain  gage  responses;  the  system  is  capable  of 
sampling  rates  up  to  50  kHz. 

The  stiffened  plate  specimen  was  subject  to  a  zero  mean  random  loading  consisting  of  a 
eomputer-simulatcd  sequence  of  10,000  endpoints  of  Rayleigh  distributed  extrema  generated 
using  a  first  order  auto-regressive  technique.  Endpoints  are  connected  by  haversine  eurves  to 
produee  a  continuous  waveform  of  unit  RMS  value.  Loading  is  customized  by  multiplying  each 
endpoint  of  the  load  sequence  by  the  product  of  desired  RMS  stress  and  average  cross  sectional 
area  of  the  specimen.  Specifically,  3  ksi  RMS  loading  was  used  in  the  tests  reported  in  this  study. 


(a)  (b) 

Figure  16:  MAHI  aluminum  hull  fatigue  specimen:  (a)  front  view;  (b)  side  view 
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Figure  17:  Strain  gage  installation  on  the  MAHI  aluminum  stiffened  hull  fatigue  specimen:  (left)  front/back 
view;  (right)  side/side  view.  Yellow  gages  are  wired  while  the  green  gages  are  wireless. 


Figure  18:  (Left)  strain  gage  installation  on  the  MAHI  specimen;  (right)  Narada  wireless  sensor  nodes  for 

data  acquisition  from  the  “green”  strain  gages. 


All  peaks  arc  followed  by  an  equal  and  opposite  trough.  All  loading  is  performed  under  load 
control.  The  10,000  point  load  sequence  was  continuously  repeated  20  to  30  times  during  day 
time  hours  and  rested  at  night.  The  stiffened  plate  specimen  was  axially  loaded  at  a  load  rate  of 
400  to  500  kips  per  second.  The  laboratory  data  acquisition  system  collected  strain  gage 
response  data  at  various  locations  on  the  aluminum  hull  specimen.  An  additional  set  of  strain 
gages  interfaced  to  the  Narada  wireless  sensor  node  is  also  installed  (many  in  locations 
redundant  to  the  location  of  the  wired  strain  gages).  Fig.  17  provides  a  drawing  of  the  test 
specimen  with  the  installed  strain  gage  sensor  layout  for  both  the  Narada  wireless  system  as  well 
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Figure  19:  Data  comparison  between  NSWCCD  tethered  system  and  Narada  wireless  system 

as  for  the  laboratory  data  acquisition  system.  Fig.  18  documents  the  test  setup.  Several  channels 
were  dedicated  to  raw  strain  gage  data  acquisition  as  well  as  the  embedded  rainflow  counting 
and  Dirlik  fatigue  life  monitoring  procedures  described  in  this  report. 

All  Narada  data  acquisition  is  conducted  at  a  sampling  rate  of  100  Hz.  Raw  strain  data 
acquisition  was  performed  continuously  as  well  as  the  embedded  rainflow  counting  (except  when 
loading  was  suspended).  Raw  strain  data  was  transmitted  from  the  Narada  to  a  receiver  base 
station  after  10,000  samples  were  eolleeted  (roughly  every  100  seconds)  and  routed  to  a  server 
(laptop)  via  a  USB  eonneetion  for  long-term  storage.  Processed  rainflow  counting  data  was 
eolleeted  periodically  for  updates.  Data  for  the  embedded  Dirlik  procedure  was  eolleeted  for  one 
minute  duration  every  hour.  Therefore,  the  resulting  stress  range  PDF  was  applied  to  a  one  hour 
time  period  when  ealeulating  the  accumulated  damage.  Processed  data  from  the  Dirlik  procedure 
was  also  transmitted  periodically.  Efforts  were  made  to  sample  the  complete  range  of  the  10,000 
point  load  cycle  so  as  to  acquire  all  aspects  of  the  load  signal. 

Results  from  strain  gages  corresponding  to  Narada  wireless  sensing  unit  channels  and  NSWCCD 
laboratory  data  acquisition  channels  are  compared  for  raw  strain  gage  response,  rainflow 
counting  histograms,  and  stress  range  PDFs.  An  effort  is  made  to  eontpare  strain  gages  in  areas 
assumed  to  experience  similar  stresses  under  uniaxial  tension-compression.  According  to  Fig.  1 7, 
one  would  expect  Narada  gage  #  8  and  NSWCCD  gage  #  8  to  have  similar  strain  gage  responses. 
Fig.  19  shows  that  the  Narada  wireless  sensing  unit  acquires  accurate  strain  data  needed  for 
reliable  fatigue  life  analysis.  After  collecting  response  data,  the  embedded  fatigue  algorithms  are 
executed.  In  Fig.  20,  histograms  for  the  rainflow  counting  show  the  cycle  extractions  of  Narada 
strain  gage  #6.  For  the  purposes  of  comparison,  this  particular  unit  transmitted  raw  strain  data  as 
well  as  the  embedded  histogram  periodically.  Fig.  20(a)  shows  the  histogram  output  from  the 
embedded  procedure  after  4  days  of  continuous  fatigue  life  processing  while  Fig.  20(b)  shows 
the  output  from  the  rainflow  counting  on  the  raw  strain  gage  data  in  post-proeessing.  Final 
damage  estimates  compare  well  but  poorly  reflect  indication  of  crack  initiation  due  to  the  use  of 
an  arbitrary  S-N  curve  in  this  study.  Future  NSWCCD  laboratory  research  for  this  project 
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Figure  20:  Rainflow  counting  histogram  on  Day  4:  (a)  embedded;  (b)  post  processing 

includes  constructing  system  specific  S-N  curves  so  that  final  damage  estimates  would  more 
aeeurately  refleet  eraek  initiation  (D  =  1). 

To  validate  the  embedded  Dirlik  procedure,  the  stress  range  PDF  produced  by  Narada  strain 
gage  #7  (calculated  and  stored  on-board)  is  obtained  and  compared  with  raw  strain  data  in  post- 
process.  Fig.  21  shows  stress  range  PDFs  generated  at  strain  gage  #7  with  the  one  hour  damage 
accumulation  from  caeh  PDF.  The  comparison  of  the  two  embedded  procedures  (Fig.  20  and 
Fig.  21)  along  with  the  aeeuraey  of  the  Narada  wireless  system  against  the  tethered  system  (Fig. 
19)  displays  the  ability  of  the  proposed  decentralized  fatigue  life  monitoring  system  to  perform 
analyses  typieally  reserved  for  eentral  processing.  Unfortunately,  damage  accumulation  is 
difficult  to  compare  as  all  strain  gage  channels  were  plaeed  at  unique  positions  on  the  specimen. 
However,  damage  accumulation  for  NSWCCD  strain  gage  #9  and  Narada  strain  gage  #9  are 
shown  in  Fig.  22  beeause  of  their  close  proximity’.  Very  similar  damage  accumulation  can  be 
observed  for  the  two  strain  gages  sensed  on  different  networks  and  processed  by  centralized  and 
decentralized  systems. 
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Figure  21:  Stress  range  PDF  obtained  by  the  Dirlik  formula:  (a)  embedded  procedure;  (b)  raw  data  in  post 

processing. 
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Figure  22:  Damage  accumulation  for  two  closely  situated  strain  gage  sensors  (NSWCCD  #9  &  Narada  #9) 
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5.  Summary  and  Future  Work 

This  report  describes  a  major  effort  towards  the  embedment  of  damage  detection  algorithms  into 
the  computational  cores  of  wireless  sensors  for  automate  data  interrogation  in  a  wireless  hull 
monitoring  system.  The  advantages  of  embedded  data  processing  include  the  preservation  of 
wireless  node  power  and  the  ability  to  deploy  larger  number  of  strain  gages  in  a  hull  without 
concern  of  data  inundation  issues.  Presented  arc  two  decentralized  fatigue  life  monitoring 
schemes  easily  embedded  and  implemented  in  a  wireless  sensor  network:  rainflow  peak  and 
mean  counting  and  the  Dirlik  spectral  methods.  The  Narada  wireless  sensing  node  was  shown  to 
be  accurate  in  its  recorded  time  history  data  when  compared  to  a  traditional  tethered  data 
acquisition  system.  Furthermore,  both  the  rainflow  and  Dirlik  procedures  were  consistent  and 
produced  fatigue  estimation  consistent  with  one  another.  The  ability  of  the  Narada  wireless 
system  to  accurately  acquire  data  and  then  accurately  process  data  provides  us  with  the  tool  for 
decentralized  processing  in  future  wireless  hull  monitoring  systems. 

Current  ship  O&M  strategies  rely  upon  frequent  visual  inspection  of  ship  hulls  to  determine 
damage.  However,  in  the  case  of  aluminum  ships,  this  task  is  more  difficult  due  to  layers  of  fire 
installation  installed  in  the  hull.  In  addition,  visual  inspections  of  the  entire  hull  are  both  costly 
and  labor-intensive.  Therefore,  the  installation  of  a  structural  health  monitoring  (SHM)  system 
coupled  with  an  effective  damage  detection  methodology  such  as  the  approach  presented  in  this 
report  can  reduce  the  cost  of  inspection  by  providing  inspectors  with  a  prioritized  list  of  probable 
areas  of  damage.  Our  future  investigation  will  extend  the  proposed  fatigue  estimation  strategics 
to  include  reliability-based  methods  of  estimating  remaining  life  in  a  hull.  In  addition,  inclusion 
of  the  fatigue  estimates  into  a  comprehensive  life  cycle  analysis  will  be  sought  to  dramatically 
reduce  the  total  cost  of  ownership  of  the  instrumented  vessel. 
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